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ABSTRACT 

Warm cores (or hot corinos) around low-mass protostellar objects show a rich chem- 
istry with strong spatial variations. This chemistry is generally attributed to the sub- 
limation of icy mantles on dust grains initiated by the warming effect of the stellar 
radiation. We have used a model of the chemistry in warm cores in which the sublima- 
tion process is based on extensive laboratory data; these data indicate that sublimation 
from mixed ices occurs in several well-defined temperature bands. We have determined 
the position of these bands for the slow warming by a solar-mass star. The resulting 
chemistry is dominated by the sublimation process and by subsequent gas-phase re- 
actions; strong spatial and temporal variations in certain molecular species are found 
to occur, and our results are, in general, consistent with observational results for the 
well-studied source IRAS 16293-2422. The model used is similar to one that describes 
the chemistry of hot cores. We infer that the chemistry of both hot cores and warm 
cores may be described by the same model (suitably adjusted for different physical 
parameters). 
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1 INTRODUCTION 

The early stages of low-mass star formation are observed 
at millimetre and submillimetre wavelengths, where both 
the dust continuum as well as molecular emission from the 
outer regions of the envelope can be detected. The dynam- 
ical aspect of low-mass star formation, with low mass stars 
forming from dense and cold condensations called prestellar 
cores, and evolving into Class and I sources, is relatively 
well defined. What is less well understood is how prestellar 
cores evolve chemically into protostellar cores which then, 
eventually, develop protoplanetary disks. Nevertheless, it is 
now clear that, as in the case of massive star formation, 
towards the end of the collapse phase the gas and dust 
are dense and cold enough that (other than H2 and He) 
material has frozen onto the dust and that surface chem- 
istry has taken place, leading to a rich gas phase chemistry 
once the star heats the surrounding dust so that icy man- 
tles returned to the gas phase. The warm cores where the 
icy mantles sublime are often called 'hot corinos' (e.g. |Cec-| 



cyanide CH3CN, and more complex species such as ethylene 
glyco l (CH20H)2 (e.g. [Cazaux et al.||2003| |Bottinelh et al. 



2004 Hollis et al. 2004 Bisschop et al. 20081. For a com- 



prehensive review on complex organic interstellar molecules 



we refer the reader to Herbst & van Dishoeck ( 2009 1. In ad- 



dition, they show some interesting deuterated species (e.g. 
[Ceccarelli et al.]|1998a[ |Caux et al.||2005i ), and long-chain 
unsaturated hydrocarbons and cyanopolyynes (Sakai et al 



2008 1. In fact, the envelope of Class sources can be divided 



carelli 2004 20071 due to their chemical resemblance with 



hot cores. Recent observations toward several warm cores 
around low-mass solar-like protostars have shown that they 
are rich in large, complex organic molecules (COMs) includ- 
ing methanol CH3OH, methyl formate HCOOCH3, methyl 
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chemically into two zones: the zone where all the ices sub- 
lime (dust temperature ^ 100 K), i.e. the warm core region, 
and the outer cold envelope (dust temperature 5C 100 K) 
(e.g. |Ceccarelli|[2005[ ). The chemistry in the outer envelope 
is similar to that in pre-stellar cores, which are cold (~ 10 
K) and dense ( ^ 1 x 10^ cm~^) condensations believed to be 
precursors of the Class sources ( [Bergin fc Tafalia||2007[ ). 
Many molecules in the outer envelope are frozen onto the 
grain mantles, so that their abundances are often lower than 
those in molecular clouds. In the inner regions, sublimation 
is clearly important. 

The first model of the thermal structure of a low- 
mass protostellar envelope was developed by |Ceccarelli et aL] 
( 1996 1. In this model, the authors considered the chemistry 



and resulting line emission within the inside-out collapse 
model ( |Shu||1977| ) but down to much smaller scales of ~ 10 
AU. In their model, the gas temperature in the innermost 
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regions of the envelope reaches ~ 100 K resulting in the 
sublimation of H2O from grain mantles. 

More recent models (e.g. Ceccarelli et al.|2000b Maret 
|et al.| |2004[ ) suggest that molecular emission (e.g. from 
H2CO) is emitted by two components: a cold H2CO-poor 
outer envelope and a warm H2C0-rich core. The model cal- 
culations show a dramatic increase in the abundances of 
H2CO, H2O, and SiO in the inner, warmer, and denser re- 
gion of the envelope. This abrupt change occurs when the 
dust temperature exceeds 100 K and grain mantles sublime. 

Evidently, in situations where a cool outer envelope and 
a warmer core exist, sublimation - and the subsequent gas- 
phase chemistry that this causes - will be important. Also, 
the sublimation will be a time- and space-dependent event 
as the central protostar warms up. Given the role of sub- 
limation, it is important to take account of recent labora- 
tory studies of the sublimation of ices ( Collings et al.|2003a 



2004 1 which show that different species enter the gas phase 



at different temperatures, and that the sublimation occurs 
not continuously through a gradual warm-up but in several 
well-defined and narrow temperature bands. Such a process 
may give rise to a rather different gas-phase chemistry than 
a process in which instantaneous sublimation of all species 
in a mantle occurs, as found to be the case for hot cores 
( Vitiet al.|2004 |. 

Our purpose in this paper is, therefore, to model the 
sublimation process, taking account of the laboratory data, 
in the context of warm cores in low-mass star-forming re- 
gions. We wish to assess whether the more complex subli- 
mation mechanisms indicated by the laboratory studies sig- 
nificantly affect the detectable chemistry in these objects, 
as compared to models in which the desorption is instanta- 
neous for all species. The model we use is developed from a 
similar model used for the study of hot cores in regions of 
massive star formation ( Viti et al.|2004 l. Further, if the two 
similar models can successfully describe cores in both types 
of star-forming region, then this may be taken as support 
for the view that star formation is carried out by essentially 
the same process, regardless of mass. We describe our model 
in Section 2, the results are given in Section 3, and we make 
some conclusions in Section 4. 



2 THE MODEL 

In order to investigate the chemical evolution of low-mass 
star forming regions, we have used a time-dependent gas- 



grain chemical model adapted from the Viti et al. (20041 
hot core model. In this work, our model follows the chemi- 
cal evolution of a free-fall collapsing cloud in the process of 
forming a low-mass star, then explores the chemical evolu- 
tion of the remnant of the warm core located in the vicin- 
ity of the newly formed star. The core is represented by 
a uniform slab subdivided into 5 shells represented by 5 
depth points of increasing visual extinction from the edge of 
the core to its centre. The results are obtained from a two- 
phase calculation: Phase I, the collapsing phase, in which 
diffuse, mainly atomic material of initial number density 
~ 400 cm""^ at a temperature of 10 K undergoes collapse 
following the so-called modified free-fall collapse described 
in [Rawlings et al.| ( |1992[ ) until the number density reaches 
the typical density of a warm core (treated as a free param- 



eter). The initial elemental abundances relative to hydrogen 
that we have adopted are listed in Table [l]( | Viti et al.|2004 |. 
During this time, gas-phase chemistry occurs and atoms and 
molecules are depleted on to grain surfaces and hydrogenate 
when possible. The depletion efficiency is determined by the 
fraction of the gas-phase material that is frozen on to the 
grains. This fraction is arranged by adjusting the grain sur- 
face area per unit volume, and assumes a sticking probability 
of unity for all species. The fraction of material on grains is 
then dependent on the product of the sticking probability 
and the amount of cross-section provided per unit volume by 
the adopted grain size distribution. Note that, as the chem- 
istry is time-dependent, different species form at different 
times, and as a consequence the material frozen out on the 
grains at any one time is representative not of the whole gas 
but of selected species. Apart from direct hydrogenation, the 
only other surface reactions we include are the formation of 
CH3OH and H2CO from CO and of CH3CN from the reac- 
tion of methane, CH4, with HCN, as it has been shown that 
gas phase reactions are not sufficient to form CII3OII and 
CH3CN (e.g. l Tielens fc Hagen|p82l [Watanabe et al.||2003 



Garrod et al.|200 8D 

In Phase II, the warming-up phase, we follow the chem- 
istry of the remnant core, after the star is born, for ~ 10^ 
yrs, assuming a uniform density throughout the core. In this 
phase, the central star heats up the surrounding gas and 
dust, causing sublimation of the icy mantles. In the present 
work, we simulate this heating effect in the same way as in 



Viti et al. ( 2004 1. We assume that the presence of an infrared 



source in the centre of the core or in its vicinity causes an 
increase in the gas and dust temperature. The temperature 
is a function of the luminosity (and therefore age) of the 

(20041, but here we adjust the 



Viti et al. 



protostar as m 
luminosity/ mass power law to account for lower protostellar 
masses than in the hot core case. The temperature reaches 
its maximum (assumed to be 100 K) at ~ 10^ years which 



is the assumed (typical) age of Class sources ( Andre et al. 
1993|. 



During the warming-up phase (Phase II) of the core, 
the mantle species desorb in various temperature bands. 
[Collings et al.| ( |2003a| |2004[ ) reported results of Tempera- 
ture Programme Desorption (TPD) experiments for a large 
number of molecules. They showed that molecular ices can 
be divided into five categories: (i) CO-like; (ii) Il20-like; 
(iii) intermediate; (iv) reactive; (v) refractory. For each cat- 
egory one can estimate the fraction of a particular molecu- 
lar species that is desorbed in various temperature bands. 
These bands are from 1) the pure species; 2) a monomolec- 
ular layer on H2O ice; 3) desorption during the amorphous- 
to-crystalline H2O ice conversion (the 'volcano' effect); 4) 
co-desorption when the H2O ice desorbs; and 5) desorption 
of the bare grain surface (see Collings et al.|[2004 for more 
details) . 

Based on the previously mentioned experimental re- 
sults, [Collings et al.| ( |2003b[ ) constructed a rate model to 
explain desorption processes in a phenomenological man- 
ner over relevant astronomical time scales. By running this 
chemical kinetic simulation of water-ice desorption, we cal- 
culated the different desorption temperatures using a power- 
law temperature profile fitted to a Sun-like star. |Viti et aL 
(|2004[ ) used the observed luminosity function of [Molinari 
et al. ( 2000 1, and correlated the effective temperature of the 
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gas with the age of the accreting protostar through a sim- 



ple power law. Nomura & Millar ( 2004 \ solved the radiative 



transfer equation for a clump warmed by a central star as- 
suming that the radial profile of the temperature is defined 
as the inverse square law T{r) oc r~^'^, which is compara- 



ble to the profile obtained by Schoier et al. ( 2002 1 in which 
T{r) ex r""*. In our model, the temperature evolves in both 
time and space. The temperature profile in our model, shown 
in Eq.(l), is derived by combining the approach of Viti et al 



1 2004 1 with that of Nomura & Millar ( 2004 1 



Td{t,d) = 10 + A{ty X [d/RY 



K 



(1) 



where Td{t, d) is the temperature profile of gas and dust 
in the core surrounding the stellar object, t is the evolution- 
ary age of the collapsing core, d is the distance from the 
core centre, and R is the core radius. A and B axe two con- 
stants derived from the boundary conditions, in this case, 
the temperature at ~ 150 AU from the star at i = (T = 
10 K) and t = lO'^ yrs (T = 100 K). The latter condition 
corresponds to the temperature obtained by [Schoier et al.| 
(20021. We find that the temperature profile of a typical 



warm core capable of reproducing with high accuracy the 
empirical assumptions we consider is given by 



T4t,d) = 10 + 0.1927(t)"-''^^'' X {d/R)-°-' 



K 



(2) 



Eq. (2), was then inserted in the chemical kinetic program 
( CoUings et al.|2003b l to estimate the position of the differ- 
ent desorption bands. From this equation it can be seen that 
the temperature rises slower than in the case of high mass 
stars ( Viti et al.|2004 \ affecting, as a consequence, the times 
at which the desorption of the different species occur. The 
desorption mechanisms, throughout the core, are time- as 
well as space-dependent processes. The closer to the heating 
source the earlier the desorption occurs. The volcano and 
co-desorptions will occur at slightly different temperatures 
(and times) depending on the depth within the core but they 
are around 84 K (at ~ 5x10* yrs, for the inner shell) and 95 
K (at ~ 7x10* yrs, for the inner shell) respectively. These 
results are in line with those previously calculated for hot 
cores where the desorption temperatures decrease as a func- 
tion of slowing down of the heating rate so that the lower 
the stellar mass, the lower the volcano and co-desorption 
temperatures (see Table 2 in [Viti et al.||200^ . 

Our chemical network is a modified version of the 



UMIST database ( |Le Teuff et aL][2000i ) with 127 gas phase 
species and 42 mantle species interacting in 1871 chemi- 
cal reactions. Photo-reactions were included, taking into ac- 
count both the external interstellar radiation field and the 
internal cosmic ray induced UV field. Both direct and in- 
direct ionisation by cosmic rays were also inclu ded, using a 
cosmic ray ionisation rate, (^, of 1.3 x 10~^^s~^ ( Lepp|l992 |. 
Observations of emission lines (e.g from H2CO and H2O) 
toward low-mass protostellar objects reveal that there is a 
gradient in both density and temperature. The dust tem- 
perature reaches 100 K in cores where the densities vary be- 
tween 1 x lO'^cm"^, such as the case of IRA S 16293-2422 (e.g 



Ceccarelli et al. 2000a 



CeccarelH 2004 1 and 2 x lO^cm" 



for NGC 1333-IRAS 4A, B ( |Maret et al.|i2004 ), and their 
radii range from 200 AU down to 27 AU. 

In order to investigate the sensitivity of the chemistry 
to the physical parameters of warm cores around low-mass 
stars, we ran a grid of chemical models where we varied (i) 



Table 1. Model initial physical parameters and elemental abun- 
dances relative to hydrogen. 



Physical parameters 



Core density f 

Core maximum temperature 

Core radius 

Depletion percentage f 



1.0 x 10'^ - 2.0 X 10* cm- 

100 K 

30 - 150 AU 

85 - 100% 



Initial elemental abundances 



Carbon 

Oxygen 

Nitrogen 

Sulphur 

Helium 

Magnesium 



1.79 X 10"* 
4.45 X 10"* 
8.52 X 10-5 
1.43 X 10"^ 
7.50 X 10-2 
5.12 X 10-*^ 



t This parameter varies only during the collapsing phase (Phase 
I). 

the final density of the collapsing core and (ii) the percentage 
of the accreted species on to grain surfaces. A characteristic 
of our chemical model is that we simultaneously follow the 
chemistry as a function of time and depth: the mantle evap- 
oration is temperature-dependent and the temperature pro- 
file is dependent on both the distance from the core centre 
and time. Table [l] summarises the various physical param- 
eters used in our model. The ice composition at the end of 
Phase I will depend both on the final density of the collaps- 
ing core and the percentage of accreted species on the dust. 
Typical abundances (with respect to the total number of hy- 
drogen nuclei) for high density and high depletion, for the 
most relevant species i.e H2O, CO, H2CO, CH4, CH3OH, are 
4.0x 10-*, 3.7x 10-^ 3.6x 10-^ 1.4x 10-*, and 3.8x IQ-'' re- 
spectively. The computed abundance of CII4 in the ice is 
high and arises as a consequence of our assumption that hy- 
drogenation (in this case, of carbon) is efficient. Detection 
by ISO of CH4 ice show lower values (van Dishocck 2004]). 
However, the ISO beam encompasses material outside the 
core and the ISO observations are of young stellar objects 
in which sublimation may have already reduced the CII4 
component. Hence, the detection may be under-estimating 
the CH4 ice abundance within the pre-stellar core. 



3 RESULTS 

3.1 Chemical trends 

Figure 1 shows how the computed fractional abundance 
(with respect to the total number of hydrogen nuclei) of 
selected species varies as a function of the radius of the core 
at around 4x 10* yrs: not surprisingly the inner parts of the 
core show an increase of the molecular fractional abundances 
for most of the species while they are very low in the outer 
regions of the core. Such trends are also present as a func- 
tion of time as shown in Figures 2 and 3, where fractional 
abundances of selected species as a function of (logarithmic) 
time at visual extinction of ~ 140 mags for a warm core at 
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two different densities, are presented (2.0 x lO^cm""^, solid 
line, and 1.0 x lO^cm"^, dashed line). 

All the species in our sample show sudden increases in 
their abundances corresponding to their multi-step sublima- 
tion from grain mantles (see Collings et al.|2004 |. However, 
the fractional abundances of H2CO and CH3OH are inter- 
esting because they show extra 'jumps' which cannot be 
explained by sublimation from grains. The increase in the 
abundances of H2CO and CH3OH around these jumps is as- 
sociated with the 'volcano' desorption of CH4 from grains. 
The increase in H2CO and CH3OH seen in conjunction with 
the CH4 'volcano' desorption is due to a sequence of gas- 
phase reactions involving CH4. Although, in our model, both 
H2CO and CH3OH are formed on grains, via the succes- 
sive hydrogenation of CO (as confirmed both experimen- 
tally and theoretically (e.g. Tielens fc Hageiil|1982 Watan- 



sulphur-bearing molecules may therefore be used as tracers 
of the physical characteristics of warm cores. The SO/CS, 
SO/H2CS, and SO/OCS ratios are less than 1 at early times 
(s; 1 X 10* yrs), while those of H2S/SO2 and H2S/H2CS are 
higher than 2 at later times (^ 6 x 10® yrs). The slowing 
down of the temperature increase allows more species to sur- 
vive in the gas-phase, and hence show higher abundances 
at late times. This result is in line with the prediction of 



both Hatchell et al. ( 1998 1 and Viti et al. ( 2004 1 for the 



change of the abundances of sulphur-bearing species. The 
key difference between low and high mass cores is the abun- 
dance of H2S as a function of time: this species decreases 
much more drastically with time for low-mass cores. Note 



that H2S drives the sulphur bearing chemistry (see Hatchell 
et al.|[l998 Viti et al.|2004l. These ratios may therefore, in 



principle, be used as evolutionary indicators for warm cores. 



abe et al!]|2003 1 there is still a contribution from gas-phase although, as Wakelam et al. ( 2004 1 point out, care needs 



chemistry, especially for H2CO. Our detailed chemical anal- 
ysis shows that at the time of the jump, the formation of 
H2CO is dominated by the oxidisation of CH3 which is in 
turn formed efficiently by the destruction of CH4 by He""" 
ions. Methanol is found to be a secondary product of CH^ 
which is mainly formed (at that time) via two ion-molecule 
reactions involving CH4; one with HJ and the other with 
N2H"'". Destruction of CH^ leads to CH|^ which reacts with 
water to form CHsOH^ which recombines to give methanol. 
These sequences of reactions which involve CH4 contribute 
to the observed amount of both H2CO and CH3OH at the 
'jump'. In fact, the 'jump' in abundance for methanol is 
rather small because the CH3OHJ recombination reaction 
is rather inefficient (see Geppert et al.|2006 1 

For most species, the abundances are lower at lower 
density with the exceptions of SO and SO2. SO is mainly 
produced by reactions involving atomic and molecular oxy- 
gen during early and late times respectively, while it is de- 
stroyed by HCO^ via radical-molecule reactions which are 
more pronounced during late times. SO is also destroyed 
while forming SO2 via oxidisation and radical-molecule re- 
actions, during both early and late times. At high densities, 
SO is destroyed more efficiently than in less dense cores. For 
less dense cores, the chemistry reaches steady state earlier 
and the evaporation peaks are more pronounced. 

In Figures 4 and 5 we show models with depletion effi- 
ciencies of 100% and 85% respectively. Lowering the deple- 
tion efficiency reduces the yield of most species. However, 
some species (Fig. 4; e.g. SO, SO2) show higher abundances 
at lower depletion efficiencies: this is due to the higher abun- 
dance of atomic sulphur in the gas phase during Phase I 
which leads to an appreciable amount of SO and SO2 even 
before the warm core forms. In general, models where all 
heavy elements are frozen out give abundances that are more 
in line with the observed abundances (see Table 2 and Sec- 
tion 3.2). We have also run a model where sulphur does not 



hydrogenate as it freezes but remains in atomic form ( Wake 
lam et aL]|2004| and, as expected, we find that this affects 



the early stages of Phase II by lowering the abundance of 
H2S by a factor of 30; however, by 10^ yrs the bulk of the 
sulphur is locked in CS, SO and SO2, regardless of its initial 
form. 

From our results so far it is clear that, as in the case of 
hot cores, sulphur-bearing species show the most variation, 
whether with time, radius, density or depletion. Ratios of 



to be taken because sulphur-bearing ratios at early times 
strongly depend on the nature of sulphur in the icy mantles. 



3.2 Comparison with hot core chemical models 

The main aim of this study was to determine whether the 
predicted chemistry of warm cores differs substantially from 
that of hot cores and in particular whether we can identify 
molecules that would be enhanced in warm cores but not 



in hot cores. Viti et al. (20041 ran chemical models of hot 



cores for different stellar masses. In this work we compare 
our results for a one solar mass warm core with those of lVitil 



et al. ( 2004 1 for hot cores. Typical hot cores (as modelled by 
Viti et al.|2004 1 have temperatures of 300 K, sizes of about 
0.03 pc, and densities of the order of 10^ cm^^; in this pa- 
per we concentrate on the inner part of a typical warm core 
(or 'hot corino') and thus employ much smaller sizes, higher 
densities and lower temperatures (see Table [l| . These dif- 
ferences in initial conditions lead to two main differences in 
the chemistry: firstly we find that in general hot cores pos- 
sess higher fractional abundances, at late times (^ 2 x 10® 
yrs), while those abundances are higher in warm cores at 
early stages. This is most likely due to the higher densi- 
ties in warn cores (leading to a richer chemistry both in 
both phases) and the different evaporation times (leading 
to different species being present in the gas phase at differ- 
ent times). Secondly, the evolutionary profiles of the species 
show almost the same trends shifted toward later times for 
hot cores. The observed jump in the abundances of species 



(such as H2CO) in warm cores (e.g. Ceccarelli et al. 1996 



|2000a|b | is not observed for the hot core case, which may 
indicate that these jumps are essential features in the warm 
core evolutionary profiles for those species. COMs show high 
abundances comparable to those found in hot cores, in par- 
ticular for CH3OH, CH3CN and C2H5OH. From the above 
discussion we conclude that the set of species used as evo- 
lutionary indicators for hot cores may be used for the same 
purpose in warm cores. 



3.3 Comparisons with observations of low^-mass 
cores 

While we do not attempt here to model any particular ob- 
served warm core, the physical parameters of one of our 
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models with density of 2.0 x 10 cm~ seem to be con- 
sistent with those observed toward the well-studied warm 
core around the solar type Class source IRAS 16392-2422 
(e.g |Ceccarelli et al.|1999 , Ccccarelli_2004) [20071 - These au- 
thors find that at a radius of 150 AU and density 2.0 x 10* 
cm~'^, the dust temperature is ~ 100 K. Hence we may 
qualitatively compare our results to the observations of 
this core. Ceccarelli et al. ( 2000a|b I recorded sudden (spa- 
tially abrupt) increases in the abundances of a few observed 
species such as H2O, SiO, and H2CO. Schoier et al. (20021 



found that in order to reproduce the physical properties 
of the dust and gas components constituting the material 
in the circumstellar envelope of IRAS 16293-2422 a jump 
in the abundance of II2CO and other species was indeed 
needed. Both of these papers proposed that these jumps, 
recorded around 80-90 K, result from the evaporation of 
mantle species in the hot-core like region around the proto- 
star. Our model with density of 2.0 x 10* cm~^ does indeed 
reproduce the jumps in the profiles of H2CO and CH3OH, 
around 84 K, confirming their findings. 

A tentative comparison between our model calcula- 
tions averaged through the core and observational fractional 
abundances in IRAS 16293-2442 is shown in Table 2, at time 
around 8.5 x 10* yrs, when the dust temperature reaches 
~ 100 K. Because our model calculations are for a core of 
radius 150 AU, we compare our results with the observa- 



tions of the inner core by Schoier et al. (2002j), apart from 
the measurement of CH3CN that was taken from |Bisschop| 
et al.|(|2008|). 



Sulphur-bearing species are, on average, in good agree- 
ment with observations during the early evolutionary stages 
(t ^ 7 X 10* yrs) and hence they are good tracers for young 
cores with temperatures Sj 84 K. On the other hand, we find 
that large complex molecules show good agreement with ob- 
servations at later times (i Js 9 x 10* yrs) where the temper- 
ature exceeds 100 K. This result supports the view of grain 
mantles as factories of complex species and it underlines 
the importance of an accurate desorption treatment in the 
chemical modelling to provide abrupt changes in molecular 
abundances. 



4 CONCLUSIONS 

We report here results from a chemical model of a warm core 
(or hot corino) around a low-mass star. The core contains 
a protostellar object of solar mass whose radiation warms 
its inner part. This leads to sublimation of icy mantles de- 
posited during the preceding collapse phase. We have intro- 
duced a detailed model of sublimation based on extensive 
laboratory results of CoUings et al. (20041. These experi- 



ments imply that sublimation from mixed ices during the 
warming process occurs in several well-defined temperature 
bands. We have computed the positions of these bands for 
the solar-mass model and follow the resulting chemistry in 
a time- and space-dependent manner. 

Although we have not attempted to model any specific 
astronomical source, the physical parameters we have used 
are sufficiently similar to those inferred from observations of 
a well-studied object, IRAS 16293-2422, for us to compare 
the model results with observations. The predicted and in- 
ferred molecular abundances agree reasonably well. 



The model that we have used is essentially a hot core 
model in which several physical parameters have been ad- 
justed to be appropriate to warming by a solar-mass pro- 
tostar. The generally satisfactory behaviour of the model 
results suggests that the chemical processes for both warm 
cores in low-mass star forming regions and hot cores in high 
mass star forming regions are similar. 

Observations of warm cores show the existence at par- 
ticular values of Av of very substantial 'jumps' in abundance 
(by several orders of magnitude) in certain molecular abun- 
dances. We find that these jumps arise as a consequence of 
sublimation and subsequent rapid gas-phase chemistry. Such 
jumps do not occur in models that describe sublimation as 
an instantaneous process, rather than the staged process we 
have introduced. 

We predict that ratios of certain sulphur species should 
be excellent tracers of physical conditions in warm cores, 
regardless of the (unknown) total sulphur abundance. 
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Figure 1. Fractional abundances with respect to the total hydrogen nuclei as a function of the core radius (Rcore) for selected organic 
species at 4x10* years for a model with density of 2x10* cm""^. 
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Figure 2. The time evolution of tlie fractional abundances of S-bearing molecules at Av = 142.7 mag. The different curves compare the 
evolution of the species at two different final densities for the collapsing cloud (see key). 
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Figure 3. The same as in Fig. (2), but for large and organic molecules. 
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Figure 4. The time evolution of tlie fractional abundances of S-bcaring molecules in a warm core region, of density nn = 2.0 X 10** 
cm~^, at a point with Ay = 142.7 mag. The profiles compare the evolution of the fractional abundances at different depletions on the 
grain surfaces (see key). 
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Figure 5. As in Figure (4) but for the large organic molecules 



